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The solar disk is a bright gamma-ray source. Surprisingly, its flux is about one order of magnitude 
higher than predicted. As a first step toward understanding the physical origin of this discrepancy, 
we perform a new analysis in 1-100 GeV using 6 years of public Fermi-LAT data. Compared to 
the previous analysis by the Fermi Collaboration, who analyzed 1.5 years of data and detected the 
solar disk in 0.1-10 GeV, we find two new and significant results: 1. In the 1-10 GeV flux (detected 
at > 5cr), we discover a significant time variation that anticorrelates with solar activity. 2. We 
detect gamma rays in 10-30 GeV at > 5a, and in 30-100 GeV at > 2a. The time variation strongly 
indicates that solar-disk gamma rays are induced by cosmic rays and that solar atmospheric magnetic 
fields play an important role. Our results provide essential clues for understanding the underlying 
gamma-ray production processes, which may allow new probes of solar atmospheric magnetic fields, 
cosmic rays in the solar system, and possible new physics. Finally, we show that the Sun is a 
promising new target for ground-based TeV gamma-ray telescopes such as HAWC and LHAASO. 

PACS numbers: 95.85.Pw, 96.50.S-, 13.85.Qk, 96.50.Vg 


I. INTRODUCTION 

The Sun is well studied and understood with a broad 
set of messengers at different energies. For example, the 
optical photon and MeV neutrino spectra confirm a de¬ 
tailed picture of the Sun as a middle-aged G-type main- 
sequence star powered by nuclear fusion [1, 2]. How¬ 
ever, the gamma-ray emission from the Sun is poorly un¬ 
derstood. Precision studies of the Sun at GeV energies 
are only now possible after the 2008 launch of the Fermi 
Gamma-Ray Space Telescope (Fermi). 

Naively, one does not expect the quiet Sun (also known 
as the steady or the quiescent Sun) to produce an ap¬ 
preciable GeV gamma-ray flux. Even though the solar 
atmospheric temperature rises to millions of Kelvin in 
the corona, it corresponds to < keV in energy. And, al¬ 
though solar flares can accelerate particles non-thermally, 
bright flares are rare and the highest-energy gamma ray 
observed from a flare is only ~ 4 GeV [3-6] . 

There are, however, two distinct processes involving 
cosmic rays that guarantee the continuous production of 
gamma rays from the vicinity of the Sun. The first con¬ 
tribution comes from the Inverse-Compton (IC) scatter¬ 
ing of cosmic-ray electrons and positrons with solar pho¬ 
tons [7-9]. The IG component appears as an extended 
halo (^ 0(10°)) around the Sun. The second contri¬ 
bution comes from the hadronic interaction of cosmic 
rays with the solar atmosphere (photosphere and chro¬ 
mosphere) [10]. The extent of this component has the 
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angular size of the Sun (~ 0.5°); we denote it (plus any 
potential non-cosmic-ray contribution) as the solar-disk 
component. 

Theoretical estimation of both components requires 
taking into account the effects of solar magnetic activity. 
Magnetic fields carried by the solar wind modulate the 
fluxes of cosmic-ray particles in the solar system [11-13]. 
This effect is expected to be stronger for the solar-disk 
component than the IC component because of the much 
closer approach to the Sun for the parent cosmic rays. In 
addition, magnetic fields in the solar atmosphere [14-16] 
affect the solar-disk component. Seckel et al. [10] (de¬ 
noted as SSG1991 in the following) showed that solar 
atmospheric magnetic fields could boost gamma-ray pro¬ 
duction through the magnetic reflection of the primary 
cosmic rays or their showers out of the Sun. Conse¬ 
quently, they estimated that the Sun could be detected 
by space-based gamma-ray telescopes. 

The first experiment to have the sensitivity to de¬ 
tect quiet Sun gamma rays was the Energetic Gamma 
Ray Experiment Telescope (EGRET) [17]. A reanaly¬ 
sis of the EGRET data later reported the first detec¬ 
tion of solar-disk gamma rays, but the flux uncertain¬ 
ties were large [18]. More recently, with the improved 
sensitivity of the Large Area Telescope (LAT) on board 
Fermi, the IC and solar-disk components were each well 
measured at 0.1-10 GeV in Abdo et al. [19] (denoted as 
Fermi2011 in the following). The IC component was de¬ 
tected out to 20° from the Sun, and was found to be 
consistent with theoretical expectations [7-9]. Although 
the observed solar-disk component satisfies the theoreti¬ 
cal upper bound derived in SSG1991 (the naive case), it 
is in complete disagreement with the nominal model of 
SSG1991, the one and only theoretical prediction: The 
observed flux is about one order of magnitude higher at all 
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energies and the spectrum shape is flatter than predicted. 
This mismatch motivates new theoretical modeling and 
new observational studies of the solar-disk gamma rays. 
The latter is the focus of this study. 

After Fermi2011, two key questions naturally surfaced 
concerning solar-disk gamma rays. First, does the solar- 
disk gamma-ray flux have a long-duration time varia¬ 
tion? In Fermi2011, after comparing to the results from 
Ref. [18], it was pointed out that a significant varia¬ 
tion of the solar-disk emission may be present. If such 
a variation is confirmed, and if it is related to the so¬ 
lar activity cycle, it could test the cosmic-ray origin of 
the gamma rays and help reveal their production mecha¬ 
nism. Second, does the Sun shine in gamma rays beyond 
lOGeV? The last two data points from the Fermi2011 
solar-disk energy spectrum suggest the spectrum might 
become softer at higher energy. Interactions of cosmic 
rays with solar magnetic fields are energy dependent; a 
spectral cutoff at high energy could reveal the end of 
magnetic field effects on the cosmic-ray interactions. It 
is only possible to answer these questions now because 
of the improved statistics and long time baseline (> 6 
years) of the Fermi-LAT data set. 

We aim to address these questions in this work, which 
is structured as follows: In Sec. II, we present our analysis 
and findings. In Sec. Ill, we first provide a short overview 
of the hadronic solar gamma-ray production by cosmic 
rays. Then we discuss future prospects for both theory 
and observation. Seasoned readers on cosmic-ray theory 
can skip the overview (Sec. Ill A) and move on to the 
rest of the section. We conclude in Sec. IV. 


II. THE SUN OBSERVED USING FERMI-LAT 
A. Outline of the Analysis 

Launched in 2008 on board Fermi, the LAT instrument 
is a pair-conversion gamma-ray detector sensitive to en¬ 
ergies from about 10“^ GeV to 10^ GeV [20, 21]. Its large 
field of view allows it to survey the whole sky. With 1.5 
years of data, Fermi2011 detected the solar-disk and IG 
components separately in 0.1-10 GeV. Since then, Fermi 
not only collected more data, but its quality has also im¬ 
proved. Fermi data are publicly available, which allows 
us to perform this study. 

Due to the apparent motion of the Sun on the sky, 
one needs to trace its position continuously with time to 
produce a Sun-centered image. Because we focus on the 
solar-disk component, all other sources of emission are 
treated as backgrounds. There are two main backgrounds 
that need to be accounted for; both are small compared 
to the signal. The first is the diffuse background that 
consists of astrophysical emission (smeared due to the 
motion of the Sun) and the detector background. The 
second background (technically, a foreground) is the IG 
component in the line of sight. Both backgrounds can be 
estimated from the data. 


We follow Fermi2011 by removing data near the Galac¬ 
tic plane and model the diffuse background using the 
fake-Sun method. In addition, we remove all the bright 
solar flares. To increase photon statistics, we relax the 
point-source cut and moon cut used Fermi2011. We 
study and take into account the possible systematics as¬ 
sociated with this step. 

To extract the solar-disk signal, we perform a likeli¬ 
hood analysis with the data binned in both energy and 
angle. This allows us to perform a simple and conser¬ 
vative analysis to characterize the main features of the 
signal. The accuracy goal of this analysis is limited by 
the systematic uncertainty of Fermi-LAT’s effective area, 
which is estimated to be about 10% [21], so we ignore 
uncertainties that are much less than that. We discuss 
possible ways to improve the analysis in Sec. III. 


B. Data Selection and Cuts 

We choose our analysis energy range to be 1-100 GeV. 
Below 1 GeV, the point spread function (PSF) of Fermi- 
LAT deteriorates rapidly, making it difficult to isolate 
the solar-disk component (in addition, the Fermi Col¬ 
laboration is performing a dedicated analysis at low en¬ 
ergies [22]). Above 100 GeV, although we find 3 pho¬ 
tons (up to ~ 300 GeV) within 1° of the center of the 
Sun in the final photon map, it is difficult to estimate 
the background contribution due to the small number of 
photons. 

We analyze the data using the Fermi science tools ver¬ 
sion v9r33p0 [23] . We use the weekly P7REP data set from 
week 010 to week 321, which covers from 2008-08-07 to 
2014-07-31. (Pass 8 data became available during the 
final stages of this work; we discuss this in Sec. III.) To 
trace the Sun’s position, we divide each week into 40 
identical time segments. Because the Sun moves ~ 7° 
per week, its positional drift per time segment is ~ 0.2°. 
This is smaller than the diameter of the Sun (~ 0.5°) 
and the LAT PSF at I GeV (~ 1°). Above 10 GeV, the 
drift becomes comparable to the PSF (~ 0.1°), which 
we mitigate by using large angular bins in the likelihood 
analysis. 

For each time segment, we adopt the standard data 
selection procedure recommended by the collaboration. 
We use gtselect to select photons from the SOURCE 
event class and to divide the events into eight energy 
bins of equal logarithmic width. We set the maximum 
zenith angle to be 100° to avoid photons coming from 
the luminous Earth limb [24, 25]. We select all photons 
within 10° of the Sun; to avoid potential edge effects, 
we define our region of interest (ROI) as a 9°-radius 
circle. The photon events are filtered using gtmktime 
with the keywords DATA_QUAL==1, LAT_C0NFIG==1, and 
ABS(R0CK_ANGLE<52). The first two keywords ensure 
that the data quality is good enough for a point-source 
analysis; the last one requires that the spacecraft be 
within the range of rocking angles used during nomi- 
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FIG. 1. Left: Stacked photon counts map of the Sun ROI in 10-100 GeV. Right: Same, but for a fake-Sun ROI (in this 
example, trailing the Sun in its path by +180 days), which is used to measure the diffuse background. The exposures of the 
two ROIs differ by < 2%. (Maps for > 0.1 GeV are shown in Fermi2011.) Visually, the solar-disk component (comparable in 
extent to the size of the Sun, as marked) is obvious; that of the IC component (decreasing with angle) is more subtle. The 
numbers of photons within 1.5° of the center are 175 versus 19; the numbers in 1.5°-9° are 844 versus 710. 


nal sky-survey observations. The filtered photon events 
are binned into photon counts maps in equatorial coor¬ 
dinates using gtbin with a pixel size 0.1° x 0.1°. The 
photon maps are stacked to construct a single map for 
each energy bin. 

To calculate the expected number of photons from 
an underlying intensity (flux per solid angle) distribu¬ 
tion, we obtain the exposure map using gt It cube and 
gtexpcube2 with identical settings as for the photon 
maps, and using the P7REP_S0URCE_V15 instrumental re¬ 
sponse function. The flux map is obtained by dividing 
the stacked photon map by the stacked exposure map. 
The total exposure in the ROI is about ~ 10^^ cm^ s, 
and is spatially uniform at the ~ 1% level in 1-100 GeV. 

To check our data selection procedures, we measure the 
gamma-ray flux from one of Fermi’s calibration sources, 
the Vela pulsar, which is the brightest steady astrophysi- 
cal gamma-ray source above 0.1 GeV. We repeat the same 
data selection procedures, except for the time segments 
used to trace the Sun, to obtain the photon map and ex¬ 
posure map. The gamma-ray flux is estimated from the 
total flux within 1.5° of Vela, after subtracting the back¬ 
ground estimated from the 6°-9° region of the same ROI. 
The flux obtained is consistent with that in Ref. [26] . 

Following Fermi2011, we remove data when jfej < 
30°, where b is the Galactic latitude. This avoids the 
bright diffuse and point-source emission from the Galac¬ 
tic plane. After this cut, the exposure time is reduced by 
~ 40% and the total photons by ~ 76%, consistent with 
the values in Fermi2011. This cut is efficient for reduc¬ 
ing background contamination, but is conservative be¬ 
cause the Galactic plane emission decreases rapidly with 
Galactic latitude. We discuss in detail the remaining 


background components in Sec. IIG. 

In Fermi20Il, data are excluded whenever a known 
point source or the Moon is within 20° of the Sun. In 
order to maximize the photon counts in high energy, we 
relax these cuts. Point sources are expected to increase 
the diffuse background by about 10%, which has mini¬ 
mal effect to our solar-disk-centric analysis. The Moon 
should not affect our analysis because its energy spec¬ 
trum falls rapidly above 1 GeV [27]. We describe in the 
next section in detail how we handle the inclusion of back¬ 
ground sources in the likelihood analysis. Imposing the 
point-source cut would reduce the exposure time by at 
least a factor of 3 (shown in Fermi2011 with IFGL), mak¬ 
ing the high-energy analysis significantly more difficult. 
(The IG component has a smaller signal-to-noise ratio. 
As a result, the point-source cut is more important for 
an IC-centric analysis, as in Fermi2011.) 

With the goal of searching for time variations in the 
solar-disk flux, we pay special attention to possible time- 
varying sources. The most important ones are solar 
flares [3-5] . During the period of bright solar flares, the 
flaring regions can emit a significant flux of gamma rays 
for a short period of time, thus contaminating the solar- 
disk signal and potentially changing the time profile of 
solar-disk flux. Only a few flares are expected to mat¬ 
ter, as solar flares are typically dim beyond a few GeV. 
Another special source is the blazar 3G 279, which over¬ 
laps the coordinates of the Sun every October [28]. This 
blazar has a flux comparable to that of the Sun and the 
Sun stays about a day near its location, hence it would 
nominally contaminate the solar-disk component at the 
~1% level. However, when it is in a flaring state, it can 
temporarily be 100 times brighter [29, 30]. We check and 
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FIG. 2. Total gamma-ray flux between 1 GeV and 1.8 GeV within 1.5° from the Sun versus time. Each bin corresponds to 
one week of observation, starting from 2008-08-07 (week 010). Periods that coincide with a bright solar flare are labeled with 
red squares; these are removed from the analyses. The horizontal grey band shows the resulting 6-year combined flux and its 
uncertainty. 


find that the Sun was never nearby during the reported 
3C 279 flares. 

Figure 1 (left) shows the stacked photon map in 10- 
100 GeV. It is clear from the density and the brightness of 
the pixels the solar disk is observed. This is the first time 
that the Sun has been detected with > 10 GeV photons. 
Gompared to the map shown in Fermi2011, which was for 
all photons above 0.1 GeV and is thus dominated by low- 
energy photons, this image is sharper due to the improved 
PSF at higher energies. The right panel of Fig. 1 is a 
“fake-Sun” photon map, used as a background estimate, 
described in the next subsection. 

Figure 2 shows the gamma-ray flux (1-1.8 GeV) within 
1.5° of the Sun as a function of time. We label the time 
periods that contain solar flares detected by Fermi at 
greater than 10 cr [31]. Some anomalously bright periods 
are correlated with solar flares, most notably the ones in 
7 March 2012 and 25 February 2014 (week 196 and 299). 
Beyond that, we do not observe any obvious excesses. 
For consistency, all labeled periods are removed from the 
Sun and fake-Sun analyses. 


C. Background Estimation 

1. Diffuse Background 

Due to the motion of the Sun on the sky, all astrophys- 
ical emission is smeared to a diffuse and isotropic back¬ 
ground. This includes truly diffuse as well as resolved and 
unresolved point-source emission. We denote this emis¬ 
sion together with the detector background (misidentihed 
cosmic rays) as the diffuse background. 

We estimate the expected contribution of the diffuse 
background in the Sun ROI using the fake-Sun method 
described in Fermi2011. We repeat identical analyses (in¬ 
cluding all cuts) at positions where the Sun would have 


been -1-60, -1-90, -1-180, and —90 days away from the actual 
time. The fake Suns traverse the same paths through the 
sky as the Sun, which allows us to measure the diffuse 
background independently. 

Figure 1 (right) shows the stacked photon map in 10- 
100 GeV for one of the fake Suns (-1-180 days). The Sun 
and fake-Sun ROIs have comparable exposures (< 2% 
difference). As a result, the small excess of photons away 
from the center of the Sun ROI already shows hints of the 
extended IC component, which becomes apparent when 
the angular distribution of the intensity is shown. 

The combination of four fake Suns allows us to esti¬ 
mate the diffuse background with better than 10% sta¬ 
tistical uncertainty. However, when comparing the in¬ 
dividual fake-Sun background estimates, we observe, at 
the low end of our energy range, ~ 10% variations among 
the fake Suns, which is larger than their individual sta¬ 
tistical uncertainties. Upon closer inspection, we found 
that this is driven by one particularly brighter fake-Sun 
ROI (-1-180), while the other three agree with each other 
at subpercent level. We check and do not find any sig¬ 
nificantly bright periods in this fake-Sun ROI. Therefore, 
this flux enhancement is likely due to one or several mild 
time-varying background sources, an arguably expected 
consequence of including point sources in the data set. 
We combine the four fake-Sun ROIs to estimate the dif¬ 
fuse background, and mitigate the potential background 
variation by adding a 10% systematic uncertainty to the 
diffuse background in the likelihood analysis. We also 
check our result using the background estimates without 
the -1-180 fake-Sun ROI. The difference is miniscule. 

We compare our combined fake-Sun background esti¬ 
mate with that from Fermi2011, and find that our back¬ 
ground estimate is higher by ^10% at the low energy 
end. Though this is consistent with systematic varia¬ 
tion described above, it could also be explained by back¬ 
ground sources. The average point-source contribution to 
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the diffuse background can be estimated using the total 
high-latitude (|6| > 20°) point-source intensity reported 
in the Fermi Isotropic Gamma-Ray Background analy¬ 
sis (see Fig. 8 in Ref. [32]). Comparing this to the diffuse 
background in the fake-Sun ROI (Fig. 3 in Fermi2011), 
point sources contribute about 10% of the total diffuse 
background, which matches the difference seen in our 
fake-Sun analysis versus that in Fermi2011. Because this 
extra small contribution affects both the Sun and fake- 
Sun ROIs, it is self-consistently modeled in the likelihood 
analysis. Nonetheless, we add an additional 10% system¬ 
atic uncertainty to the diffuse background in the likeli¬ 
hood analysis. These systematic uncertainties (20% of 
our fake-Sun estimate) take into account all the poten¬ 
tial systematics introduced to the diffuse background by 
including the point sources. 

Lastly, the gamma-ray intensity of the fake-Sun ROIs 
are found to the uniform in radial direction. This is con¬ 
sistent with the finding from Fermi20Il, which showed 
that the only source of anisotropy is the Galactic plane, 
which we have removed. This angular dependence al¬ 
lows us to separate the diffuse background from the signal 
and the IC component. 


2. Inverse Compton Emission 

In addition to the diffuse background, the extended 
IG component also contributes to the total emission in 
the Sun ROI. We model the IG component background 
using its distinctive angular distribution. Assuming the 
cosmic-ray electron density is homogeneous throughout 
the solar system, the IC component intensity is simply 
proportional to the column density of solar optical pho¬ 
tons [7-9]. This description was found to be reasonable 
in Fermi2011, especially for gamma-ray energies above 
1 GeV. With this assumption, we can approximate the 
IC intensity as cx where a is the angular distance 
from the Sun. This distribution deviates from the true 
one [7-9] slightly at large angles, and is accurate at the 
^ 5% level at the edge of our ROI. In the angular re¬ 
gion of the solar disk, the IC component is suppressed, 
which we take into account in the analysis (described be¬ 
low). Overall, small uncertainties of the shape of the IC 
component do not affect our results, as it is subdominant 
compared to the solar-disk emission in the inner 1.5°. 


D. Solar-Disk Flux Spectrum 

We use a multicomponent fit to extract the solar-disk 
component. This exploits the facts that the Sun is spa¬ 
tially concentrated (see Sec. HID for discussion on re¬ 
solving the Sun), the IC component is extended with a 
characteristic profile, and the diffuse background is spa¬ 
tially uniform. The angular information allows us to fit 
the components individually for each energy bin, without 
requiring any assumptions about the energy spectra. 


We divide the Sun ROI into angular bins that are con¬ 
centric rings of 1.5° width. This choice is guided by the 
PSF of Fermi-LAT, which is 0.8° at 1 GeV (68%). Be¬ 
cause the PSF improves above 1 GeV and flattens out by 
~ 10 GeV, the 1.5° bin ensures that the solar-disk com¬ 
ponent is always fully contained in the first angular bin. 
This criterion significantly simplifies the analysis. More¬ 
over, our choice of the uniform 1.5° angular bin across 
all energies is conservative. The PSF of Fermi-LAT im¬ 
proves at high energies, so in principle one can afford a 
smaller angular bin at higher energy bins. However, we 
expect the improvement from such an analysis will be 
small, given that the diffuse background is small. For 
simplicity, we use constant angular bins across all energy 
range. 

With this angular binning, the distribution of the 
gamma-ray flux in the Sun ROI is modeled independently 
for each energy bin, as follows: 

Si = si Sii (1) 

bf = 

j 

3 

where Si, and are the modeled photon counts 

for the solar-disk signal, as well as the IC and diffuse 
backgrounds in angular bin i. Sij is the exposure for 
a pixel j in bin i (with unit [cm^ssr]), and aij is the 
angular distance from the center to a pixel j in bin i. 
The solar-disk component is described by a Kronecker 
delta function, Sn, which indicates that the first angular 
bin fully contains the solar-disk flux. The IC compo¬ 
nent is described by a normalization factor, times 
the total exposure weighted by q~^. At the region of 
the solar disk (a < 0.27°), the IC component is strongly 
suppressed due to the anisotropy of the solar radiation 
and the occultation of the Sun [7-9]; we set the IC com¬ 
ponent to be zero in this region accordingly. The diffuse 
background component is radially isotropic, so it is only 
a normalization factor, weighted by the total ex¬ 

posure. 

For each energy bin, we perform a profile likelihood 
analysis [33, 34]. The likelihood function is a function of 
the signal parameter, si, and the nuisance parameters, 
/IC and /iiKC; 

C(K,C. /™°) = G(/™°)n 

i 

( 2 ) 

where P is the Poisson probability for the model to yield 
the observed number of photons, di. The product is taken 
over all angular bins. The Gaussian term, con¬ 

strains the diffuse background from deviating too much 
from the value determined from the fake-Sun method. 
We take the variance of the Gaussian to be 20% of the 
combined fake-Sun flux estimate, and assume that it is 
uncorrelated between energy bins. The 20% systematic 
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FIG. 3. Left: Angular distribution of the integrated intensity from 1-10 GeV in the Sun ROI. Black points show the observed 
data with statistical uncertainties only. Golored histograms show the fitted results for the signal and two backgrounds (the 
estimate of the diffuse background incorporates independent data from the fake-Sun ROIs). The inset shows the same data 
with smaller angular bins, but without the two solar components (note the different vertical scale). Right: Same, but for 
10-100 GeV (note the lower flux). 


uncertainty conservatively combines the 10% variations 
among the individual fake Suns and the 10% difference 
we observe from our fake-Sun method compared to that 
from Fermi2011. The best-fit diffuse background normal¬ 
ization in the Sun ROI is found to be within 10% of our 
fake-Sun estimate for all energy bins, which shows that 
the fake-Sun estimate is accurate and the choice of 20% 
variance for G(/®^*^) is conservative. The normalization 
of the IC component is conservatively set as a nuisance 
parameter. The final uncertainty of the extracted solar- 
disk component therefore includes the maximum normal¬ 
ization uncertainty of the IC component. 

Figure 3 shows the angular distribution of the intensity 
in coarse energy bands, given by the number of photons in 
each angular bin divided by the total exposure. The data 
points represent the total observed intensity with statisti¬ 
cal error bars only, and the colored histograms represent 
the fit for the three individual components. This simple 
model describes all features of the data well, and it is 
evident that the solar-disk component has a high signal- 
to-noise ratio. 

For each energy bin, we obtain the best-fit model pa¬ 
rameters by maximizing the likelihood function with re¬ 
spect to all model parameters. The uncertainty of the 
extracted solar-disk signal is found using the profile like¬ 
lihood function, which is the likelihood function maxi¬ 
mized over only the nuisance parameters. Assuming the 
signal parameter is Gaussian-distributed, the l-cr error 
bar of the signal is determined by where the log-profile 


likelihood function differs from the best-fit value by 1/2. 
This uncertainty determination procedure is exact when 
the sample size is large, but is found to be reasonable for 
fairly small sample sizes [34]. We check explicitly that 
the log-profile likelihood function behaves close to the 
expected parabolic shape, which verifies the Gaussian- 
distribution assumption. In addition to the uncertain¬ 
ties estimated above, the gamma-ray flux has an overall 
10% systematic uncertainty from the effective area of the 
Fermi-LAT. 

We check our result using the same 1.5-year time pe¬ 
riod as in Fermi2011. We find that our solar-disk compo¬ 
nent is consistent with that of Fermi2011, despite using 
different data sets (Pass 6 vs Pass 7), different energy 
and angular binning, different cuts, and a different anal¬ 
ysis method. This supports our analysis choices. 

For the full 6-year data set, we obtain a non-zero 
solar-disk signal in all eight energy bins from the like¬ 
lihood analysis. The detection significance can be esti¬ 
mated from the test statistic (TS = 2Alog£), given by 
two times the difference between the best fit log-profile 
likelihood function and the one with the null hypothesis 
(si = 0). The Gaussian significance, to good approxima¬ 
tion, is given by VTS [34]. As a cross check, we obtain 
comparable best-fit parameters and uncertainties using a 
simple and analysis. 

Table I summarizes our results, listing the energy bins, 
the total photon counts, and the best-fit numbers of pho¬ 
tons in the solar-disk component, and v^TS. We find that 
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FIG. 4. Energy spectrum of the solar-disk flux. Blue squares 
and statistical uncertainties (systematic uncertainties, not 
shown, are ~ 10%) are the results of our analysis with 6 
years of data. Black dots and combined statistical and sys¬ 
tematic uncertainties are the Fermi2011 results with 1.5 years 
of data. The green band shows the predicted flux range from 
the SSG1991 nominal model. 


TABLE I. Eor each energy bin, as defined, the total number 
of photons within 1.5° of the center of the Sun, the rounded 
best-fit number of photons due to the solar-disk signal, and 
the significance (VTS) of the solar-disk flux detection. 


Energies [GeV] 

Total cts. 

Best-flt solar-disk cts. 

^TS 

1.0-1.8 

1468 

961 

20.5 

1.8-3.2 

914 

628 

17.7 

3.2-5.6 

448 

329 

13.6 

5.6-10 

188 

133 

8.5 

10-18 

92 

67 

6.7 

18-32 

55 

42 

5.9 

32-56 

16 

10 

2.6 

56-100 

12 

7 

2.3 


the solar-disk component is significantly detected (> 5cr) 
up to ~ SOGeV, and is detected (> 2cr) in each of the 
last two energy bins that go up to 100 GeV. The lower 
detection significance at > SOGeV is mainly due to not 
having enough statistics to distinguish the IC and solar- 
disk components. We discuss the total solar gamma-ray 
flux more in Sec. Ill E. 

Figure 4 shows the energy spectrum of the solar-disk 
component obtained in our 6-year analysis with 1-a error 
bars. The spectrum extends to 100 GeV without an ob¬ 
vious spectrum cutoff, though for energies > 30 GeV, the 



FIG. 5. Energy spectrum of the solar-disk flux, separated into 
three periods, each of two years. The solar disk flux from first 
two years is consistent with Fermi2011, while the 1-10 GeV 
data shows a significant reduction in later periods. 


error bars are large. The spectrum can be roughly de¬ 
scribed as a single power law, oc though the power- 

law fit is not particularly good. For comparison, we also 
show the solar-disk component found in Fermi2011 and 
the SSG1991 nominal model prediction on the solar-disk 
flux, where in the former the error bars include system- 
atics. Comparing our result to that of Fermi2011, our 
analysis yields a similar spectrum with a lower normal¬ 
ization in the overlapping energy range. We find that 
this is because the underlying flux has a significant time 
variation, as detailed in the next subsection. 

Compared to the central value of the SSG1991 pre¬ 
diction, our 1-10 GeV result is still higher by a factor of 
about 5. The flux normalization of the solar-disk gamma- 
ray flux remains am unsolved puzzle. To provide more 
context on the physical implications of this disagreement, 
we discuss and provide more details about the SSG1991 
model in Sec. Ill A 2. 


E. Time Variation of the Solar-Disk Flux 

Figure 5 shows the solar-disk gamma-ray flux energy 
spectrum obtained from our analysis when we divide the 
whole data set into two-year segments (52 weeks per 
“year”). In 1-10 GeV, a decreasing trend in flux is clearly 
observed. Above 10 GeV, the situation is unclear, due to 
the large error bars. The time modulation of flux above 
GeV is already hinted at in Fig. 2, where the 1-1.8 GeV 
data showed a slow decline over the course of 6 years. 



















































FIG. 6. Left: For several sources, the ratio of the 1-10 GeV flux in each year to its 6-year average (a time-independent 
source would fluctuate around unity). The solar-disk component (blue squares) demonstrates a clear decreasing trend and 
anticorrelation with the smoothed sunspot number, a tracer of solar activity. Other sources (points displaced for clarity) should 
be and are consistent with being time-independent; see the text for details. Right: Same, but for 10-100 GeV. No obvious 
trend is observed for the disk component, but the uncertainties are large. 


To better quantify the time variation observed in 1- 
10 GeV, we first combine the data into two broad energy 
bands (1-10 and 10-100 GeV), and then find the flux 
ratio for each energy band, which is the integrated flux 
in each year relative to that averaged over 6 years. A 
time-independent source would fluctuate around unity. 

Figure 6 shows the flux ratios in these two energy 
bands. In the 1-10 GeV band, the solar-disk component 
demonstrates a significant time variation, an overall de¬ 
creasing trend, in which the extremes differ by about a 
factor of 2 to 3. We estimate the statistical significance 
of the time variation by testing the data against the null 
hypothesis (the underlying distribution is time indepen¬ 
dent) using a simple y-square test. The y-squares are 
104 and 1.6 for 5 degrees of freedom for the 1-10 GeV 10- 
100 GeV band, respectively. This shows that the time- 
variation in the 1-10 GeV band is highly significant, while 
the 10-100 GeV data is consistent with being time inde¬ 
pendent. 

We note that in Fig. 5, it can be seen that our 1-2 
year result is compatible with the Fermi2011 spectrum 
in the overlapping energy range. Given that flux only 
slightly decreases from the first year to the second year, 
this shows that our analysis with 18-months of data is 
compatible with that of Fermi2011. 

To make sure that the observed time variation is physi¬ 
cal, we check the flux ratios of several gamma-ray sources 
as control samples. First, we consider one of the fake 
Suns (-1-180 days). We find the total gamma-ray flux 
within 1.5° of its center, as in our solar-disk analysis. 


This allows us to investigate possible fluctuations of the 
diffuse background. For both energy bands, we find that 
they are consistent with being time independent. Similar 
results are obtained when other fake Suns are used. 

Second, we consider the gamma-ray flux from the Vela 
pulsar (a constant gamma-ray source), which we use to 
validate our data selection procedure in Sec. IIB. This 
allows us to check for unknown systematics in data se¬ 
lection. The flux ratios of Vela demonstrate very small 
deviations from unity in both energy bands. 

Third, we consider the total flux in the 3°-9° region 
from the Sun ROI, which allows us to check for peculiar¬ 
ities in the Sun ROI. The flux ratios are again consistent 
with being time independent for both energy bands. 

None of the control samples demonstrates any system¬ 
atic effects. This means that the observed signal time 
variation is robust, and is a feature of the underlying 
gamma-ray production processes. This variation and its 
amplitude was never quantitatively predicted and this is 
the first time it is clearly observed. 

F. Anticorrelation of the Solar-Disk Flux with 
Solar Activity 

We check whether the observed time variation is re¬ 
lated to solar activity. Our analysis period coincides 
with solar cycle 24, which started with the solar mini¬ 
mum in 2009 and reached the solar maximum in 2014. 
In Fig. 6 we overlay the yearly smoothed sunspot num- 




































9 


ber [35], which is a tracer of solar activity. Though the 
sunspot number and the solar-disk gamma-ray flux vary 
with different amplitudes, the trends are clearly opposite. 
In other words, the solar-disk gamma-ray flux anticorre¬ 
lates with solar activity at least during the first half of 
the solar cycle 24. 

This trend is also qualitatively consistent with the 
EGRET observation. The flux measured by Ref. [18] 
used data collected during 1991-1995, which is approxi¬ 
mately the second half of solar cycle 22, when solar ac¬ 
tivity was declining from the solar maximum. The anti¬ 
correlation explains the smaller flux observed by Ref. [18] 
compared to Fermi2011, who used data mainly from the 
solar minimum. 

Before this work, there was no direct evidence showing 
that the solar-disk gamma-rays are of cosmic-ray origin. 
Though only rare solar flares are found to accelerate par¬ 
ticles beyond 1 GeV, it may be possible that some yet- 
unknown solar processes continuously accelerate particles 
up to the multi-GeV energy range. However, one expects 
these solar processes would be correlated with solar ac¬ 
tivity, the opposite of the cosmic-ray framework (detailed 
in the next section). The anticorrelation with solar ac¬ 
tivity found in the solar-disk gamma-ray flux therefore 
strongly indicates that the bulk of the gamma-ray flux is 
induced by cosmic rays. (Exploration of theoretical pos¬ 
sibilities for the Sun itself to generate gamma rays that 
mimic the observed time variation is beyond the scope of 
this work.) 

It is interesting to put the amplitude of this time vari¬ 
ation into perspective, assuming the cosmic-ray produc¬ 
tion mechanism. The progenitors of 1-10 GeV solar-disk 
gamma rays are ~ 10-100 GeV cosmic-ray protons. The 
time variation (or modulation) of the cosmic-ray flux at 
Earth is known to anticorrelate with solar activity, in 
the same sense as the solar-disk gamma-ray flux found in 
this work. The cosmic-ray flux modulation at Earth is 
frequently described by the force field model with a sin¬ 
gle empirical parameter, the force field potential 4) [36- 
38]. The value of $ can be extracted from precision 
ground-based neutron observations [39-41]. We obtain 
the corresponding values for our observation period by 
averaging over the monthly values. In our analysis pe¬ 
riod, the maximum value of $ was 630 MV in 2014 and 
the minimum was 300 MV in 2009. Taking these values, 
the maximum cosmic-ray flux is larger than the mini¬ 
mum by about 15% at 10 GeV and 2% at 100 GeV. (For 
comparison, the extreme yearly values from 1964 to 2014 
are about 1200 MV and 270 MV, which corresponds to 
about 50% and 5% difference in the cosmic-ray flux am¬ 
plitude). This is too small to explain the amplitude seen 
in Fig. 6. This suggests that one needs additional mod¬ 
ulation of the cosmic-ray flux in the inner solar system, 
variations in solar atmospheric magnetic fields that can 
affect cosmic rays of such high energies, or perhaps both 
to explain the observed variation amplitude. 

In fact, the Tibet air shower array found time variation 
in observations of ^ 10 TeV cosmic-ray shadows of the 


Sun. During the solar maximum, the cosmic-ray shadows 
are shallower than during the solar minimum [42]. This 
can be explained by coronal magnetic fields: cosmic rays 
are more severely deflected by the solar atmospheric mag¬ 
netic fields during the solar maximum [42]. This implies 
that it is more difficult for cosmic rays to go deep into 
the solar atmosphere during the solar maximum, which 
is consistent with our solar-disk gamma ray observations. 

The observation of time variation in the solar-disk 
gamma-ray flux therefore provides strong support for the 
cosmic-ray framework, which we discuss in detail in the 
next section. 


III. THEORETICAL OVERVIEW AND 
OBSERVATIONAL OUTLOOK 

In this section, we first review the cosmic-ray frame¬ 
work, i.e., how solar-disk gamma rays can be produced 
from cosmic-ray interactions with the solar atmosphere. 
Experienced readers can skip the first part, and move on 
to the bulk of this section, where we discuss some future 
prospects on solar gamma-ray theory and observations. 


A. Physics of Solar-Disk Gamma Rays—The 
Cosmic-Ray Framework 

1. General Considerations 

The physics involved in the production of solar-disk 
gamma rays is complicated by the effects of magnetic 
fields. To gain some physical insights, we describe some 
general cases and approximations, following SSG1991. 

Cosmic-ray propagation from the interstellar medium 
to the surface of the Sun is known to be affected by solar 
magnetic fields carried by the solar wind. As a result, 
this propagation is also affected by solar activity [36, 37]. 
Generally, cosmic rays with energy <10 GeV observed 
at the Earth are more suppressed when the Sun is more 
active. Additional modulation of cosmic rays may occur 
when they propagate from the Earth to the Sun. 

Once cosmic rays reach the Sun, their motion is dom¬ 
inated by the magnetic fields in the corona and photo¬ 
sphere. The Larmor radius of cosmic rays near the sur¬ 
face of the Sun sets a reference energy scale, Ec- For 
cosmic-ray protons, taking the typical solar magnetic 
field strength, B ^IG, and setting the Larmor radius, 
L, to be the solar radius, Rq ~ 7 x 10® km, yields 

B, = 2xlO‘Gev(T)(^). (3) 

A similar scale is obtained for sunspots, where the length 
scale is about 10® times smaller, but the field strength is 
roughly 10® times stronger. The range of Ec was esti¬ 
mated in SSG1991 to be between ~ 3 x 10^ GeV and 
~ 2 X 10"* GeV. This scale separates the physics into 
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three regimes: Ep^ Ec, Ep E^, and Ep ~ Ec, where 
Ep is the primary cosmic-ray energy. 

When Ep 3> Ec, one can ignore the magnetic fields. 
Cosmic rays and their interaction products travel in 
straight trajectories following the initial cosmic-ray mo¬ 
mentum. In this case, only gamma rays from the Sun 
limb are observable. The Sun limb is the thin layer of 
the outer solar atmosphere that has high enough col¬ 
umn density for cosmic rays to interact, but not so much 
that gamma rays are unable to escape. This corre¬ 
sponds to a column density of 0{1) hadronic interaction 
length, which is similar to the photon absorption length. 
The Sun limb component is non-zero but is argued in 
SSG1991 to be small; it should also inherit the primary 
cosmic rays’ spectral index (~ 2.7). The Sun limb com¬ 
ponent is expected to be time-independent. 

When Ep Ec, cosmic rays propagate along solar at¬ 
mospheric magnetic field lines. Inward-pointing (towards 
the Sun) cosmic rays are funneled into magnetic flux con¬ 
centrations (or flux tubes) in the photosphere, where the 
field strength is stronger and the matter density is higher. 
Assuming adiabatic invariance, the inward-moving cos¬ 
mic rays would be reflected by the magnetic field strength 
gradient (magnetic reflection). It is then possible for the 
cosmic rays to interact with the solar atmosphere on their 
way out and to produce gamma rays that point toward 
Earth. This mechanism, suggested in SSG1991, allows 
the whole solar disk to be involved in gamma-ray pro¬ 
duction, and thus enhances the flux. Because the effects 
of magnetic fields on cosmic-ray propagation are energy 
dependent, the spectral index of the resultant gamma-ray 
flux could deviate significantly from that of the primary 
cosmic-ray spectrum. During solar maxima, the strength 
of solar atmospheric magnetic fields increases [43] , so the 
magnetic reflection of cosmic rays are expected to occur 
at higher altitudes, where the density is lower. This de¬ 
creases the gamma-ray production efficiency during solar 
maxima compared to that during solar minima, which is 
qualitatively consistent with the time variation observed 
in this work. 

When Ep ^ Ec, no simple approximation can describe 
the physics. The corresponding gamma-ray energy at 
~ O.lEc marks the transition from the low-energy regime 
to the high-energy regime. In other words, the gamma- 
ray flux, spectral index, and time-dependence should be 
intermediate between those of the two regimes above. It 
is interesting to note that the robust detection of the 
solar-disk component at 30 GeV and the non-observation 
of a spectral break in this work already requires that Ec ^ 
300 GeV, which is close to the lower bound estimated by 
SSG1991. Interestingly, the result from the Tibet air 
shower array shows that cosmic rays at ~ 10 TeV are 
still affected by solar atmospheric magnetic fields [42]. 


2. The SSG1991 Model 

We now briefly describe the SSG1991 “naive” and 
“nominal” cases for the solar-disk gamma-ray flux (see 
Ref. [10] for details). The SSG naive calculation ignores 
all the propagation and magnetic-field effects, assumes 
100% efficiency for cosmic-ray absorption in the solar 
surface, and counts all the gamma rays produced. The 
naive case, therefore, is an robust theoretical upper limit 
on how much solar-disk gamma rays can be produced by 
cosmic rays. It is not a physical model, and hence, it is 
not surprising that our flux and that from Fermi2011 is 
lower than this bound. 

The appropriate comparison with data is using the 
SSG nominal model, shown in Fig. 4. In this case the 
cosmic-ray propagation was treated as a diffusion prob¬ 
lem from the Earth to the Sun. Primarily concerning 
the Ep <C Ec case, all cosmic rays were assumed to land 
on magnetic flux tubes, and then reflected with some 
efficiency. With a chosen set of diffusion parameters, 
the cosmic-ray absorption rate was determined, which 
is roughly 0.5%. Finally, the magnetically enhanced 
gamma-ray flux was obtained by integrating the gamma- 
ray yield with the absorption rate and the path length 
distribution. The upper edge of the green band in Fig. 4 
corresponds to the extreme case where all the cascade 
products are charged and contribute to the gamma-ray 
production. The lower edge corresponds to the conser¬ 
vative case where all the cascade products are neutral, 
hence only primaries that interact after being reflected 
can contribute to the gamma-ray flux. These two cases 
bracket the theoretical uncertainty concerning the cas¬ 
cade development inside the flux tubes, but not other 
model ingredients. 

B. Prospects for Solar-Disk Gamma-Ray Theory 

As already discussed in Fermi201I, the SSG nominal 
model is unable to explain the observed gamma ray data. 
Our result, even if taken at solar maximum, is still in¬ 
consistent with the SSG nominal model. Therefore, it is 
necessary to revisit the modeling of the comic-ray frame¬ 
work. Most likely, new implementations of cosmic-ray 
physics and solar physics are needed. We will provide 
new theoretical investigations in our forthcoming papers. 

There are several key observations that the new model 
needs to address. First, it needs to reexamine the effec¬ 
tiveness of solar magnetic fields in enhancing the gamma- 
ray flux at Ep ^ Ec- In particular, SSG199I estimated 
~ 0.5% of the total available cosmic-ray energy at the 
Sun is converted to gamma rays, but observations sug¬ 
gest ~ 5%, modulo the time variation. Second, the high- 
energy gamma rays found in this work demand a proper 
treatment of the Ep ^ Ec and Ep ^ Ec regimes. Third, 
the time variation found in this work, as well as that from 
the Tibet air shower array, show that the model should 
track the variations of solar magnetic activity. Lastly, 
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the model needs to quantitatively explain the observed 
amplitude of the time variation. 

With an accurate model of gamma-ray production, 
solar gamma-ray observations can be used to constrain 
model ingredients and parameters, thus providing a new 
probe of solar atmospheric magnetic fields and of cosmic- 
ray propagation in the solar system. This is particu¬ 
larly promising given that many current and future in¬ 
struments will have excellent sensitivity for continuously 
monitoring solar gamma rays. 

With a sufficient understanding of the solar-disk 
gamma rays, it will be possible to use the Sun as a lab¬ 
oratory to test new physics. For example, a popular 
dark matter candidate is the Weakly Interacting Mas¬ 
sive Particle (WIMP), which can accumulate and anni¬ 
hilate in the core of the Sun after being gravitationally 
captured ( [44-46] ; see Ref. [47] for a recent review). Typ¬ 
ical WIMPs captured in the Sun generate negligible elec¬ 
tromagnetic signals [48]. However, non-minimal physics, 
such as inelastic dark matter [49-51] and metastable me¬ 
diators in the dark sector [52-55], can significantly en¬ 
hance the electromagnetic signatures [56-59]. Under¬ 
standing the standard model predictions is necessary to 
uncover or interpret any potential signatures from dark 
matter [10, 60-64]. For example, both the spectral infor¬ 
mation and time variation can be useful model differen¬ 
tiators. We will further discuss the implications of high- 
energy solar observations for new physics in our forth¬ 
coming papers. 


C. Prospects for the Inverse Compton Component 

In our analysis, the IC component is treated as a back¬ 
ground. However, with new data releases, which improve 
both statistics and data quality, a more precise study of 
the IC halo component is also warranted. A minor ten¬ 
sion between the data and the prediction for the IC com¬ 
ponent was found in Fermi2011, where the data seemed 
to be higher at lOGeV than expected. A more precise 
measurement is needed to clarify the situation. 

A new study of the IC component will allow one to 
use gamma rays to probe the cosmic-ray electron density 
in the solar system [8]. This is because the IC intensity 
is the product of the electron density and the photon 
density along the line of sight, with the latter being a 
known quantity. The IC component is therefore sensi¬ 
tive to electron densities from fairly close to the Sun to 
beyond the Earth’s orbit. In addition, if there is time 
variation in the IC component, its broad angular distri¬ 
bution may allow one to test the variation amplitude as a 
function of the distance to the Sun. These observations 
can help with understanding cosmic-ray modulation in 
the solar system, which despite many years of effort, is 
still under active investigation [11, 38, 65]. This approach 
is complementary to solar-disk gamma-ray observations, 
which are strongly affected by the conditions of the solar 
atmosphere. 


Similar to our analyses, it is also interesting to charac¬ 
terize the IC component beyond 10 GeV as well as search 
for long-term time variations. Because point sources are 
not removed, our analysis is not optimized for the IC 
component. With this caveat, we check the best-fit IC 
amplitude from our analysis, and we find no obvious time 
variation (only ~ 20% scatter around the mean). A more 
careful analysis is needed to provide a definitive state¬ 
ment. Analyzing the IC component is difficult at high 
energies, where statistics are low, and equally challeng¬ 
ing at low energies for Fermi-LAT, where the PSF is ^ 
10° at 100 MeV. 


D. Prospects for Fermi and Future Space Missions 

In this work, we use a straightforward analysis to char¬ 
acterize and robustly detect important features of the 
solar-disk gamma rays. Future analyses and observa¬ 
tions, with more optimized analysis procedures and im¬ 
proved data sets, can yield more precise measurements 
or even discover new features. Below, we discuss some 
possible analysis improvements with Fermi. 

At high energies, where statistics are low, one can use 
an unbinned analysis to fully utilize the information car¬ 
ried by each photon. In particular, better angular reso¬ 
lution at high energies may allow one to resolve the solar 
disk and locate hot spots (as for solar flares [4, 5]). On the 
other hand, the improved angular resolution also means 
that the solar disk can no longer be treated as a point 
source. One needs to take into account the fact that the 
astrophysical diffuse background and the IC component 
are reduced toward the solar disk [9]. This also means 
that the one should avoid using the stacking procedures 
performed in this work, which slightly smears the posi¬ 
tion of the Sun according to the length of each time seg¬ 
ment. Instead, one should select the events and calculate 
the exposure in a solar-centric coordinate system. 

For improving statistics, one can potentially develop 
more optimized cuts. For example, it is likely that the 
Galactic plane cut employed in this work can be im¬ 
proved, given that the Galactic plane gamma-ray inten¬ 
sity drops rapidly with latitude and can in principle be 
modeled. This may improve the statistics by about a fac¬ 
tor of 2. In addition, the new Fermi data release. Pass 
8 [66] , has a larger effective area and better angular reso¬ 
lution. Improving the statistics is particularly important 
for high-energy observations. 

Next-generation space gamma-ray telescopes can fur¬ 
ther improve the solar-disk observations in both time 
and energy range. The apparent anticorrelation between 
the solar-disk gamma-ray flux and solar activity suggests 
that the flux should start to increase as we start to leave 
the solar maximum. This can be checked with near¬ 
future data from Fermi. Next-generation instruments, 
such as DAMPE [67], GAMMA-400 [68], and HERD [69], 
will allow the Sun to be monitored at the GeV range even 
beyond Fermi’s lifetime. Though in principle Fermi is 
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sensitive to gamma rays down to 10 MeV, extracting the 
solar-disk signal is difScult due to the broad PSF. Future 
missions such as PANGU [70] and ComPair [71] can pro¬ 
vide improved sensitivity in the MeV range. Low-energy 
observations could provide additional information on the 
time variation and probe potential leptonic components 
or even new solar-disk gamma-ray emission mechanisms. 

E. Prospects for Ground-Based Telescopes 

To expand solar gamma-ray observations into the TeV 
range and beyond, large ground-based experiments are 
required. It is impossible for air-Cherenkov telescopes 
to observe the Sun due the bright optical emission from 
the Sun itself. The Sun, therefore, is a unique tar¬ 
get for water-Cherenkov telescopes such as HAWC and 
LHAASO. 

To assess whether water-Cherenkov telescopes can de¬ 
tect the Sun, we consider the total solar gamma-ray flux, 
including both the solar-disk and IC components. We 
estimate this flux by finding the total flux within 1.5° 
of the Sun and subtracting the diffuse background. In 
this case, the Sun is detected at > 5 cr in all eight en¬ 
ergy bins. Assuming a single power-law spectrum, the 
total solar gamma-ray flux can roughly be described by 
3.5 X 10-®(A/GeV)-2-3GeV-icm-2s-i in 1-100GeV. 

Figure 7 shows the total solar gamma-ray flux, the 
solar-disk-only component from Fermi2011, the solar- 
disk-only component found in this work, and the dif¬ 
fuse background within 1.5° of the Sun. The total solar 
gamma-ray flux is clearly much larger than the diffuse 
background. For comparison, we show also the sensitiv¬ 
ity of HAWC [72] and LHAASO [73, 74]. If the total 
solar gamma-ray flux follows the same spectral index to 
the TeV range, both HAWC and LHAASO should be 
able to detect the Sun. 

The water-Cherenkov telescopes are in a unique posi¬ 
tion to probe solar gamma rays. In particular, they are 
sensitive to the Ep ^ Ec and Ep 3> Ec regimes. Either 
a detection or an upper limit from the water-Cherenkov 
telescopes can provide valuable information on gamma- 
ray production from the Sun. 


IV. CONCLUSIONS 

Despite being the nearest star to us, much about the 
Sun’s gamma-ray emission is still poorly understood. 
Previous study by the Fermi collaboration, who used 1.5 
years of data, precisely detected the solar-disk gamma 
rays in 0.1-10 GeV. However, the flux is about ten times 
brighter than predicted. Motivated by this puzzle, we 
focus on the solar-disk component, and use 6 years of 
public Fermi data to gain a better understanding of these 
gamma rays. We employ a straightforward and conserva¬ 
tive analysis to search for new features in the gamma-ray 
flux. 



FIG. 7. Energy spectrum of gamma rays from the Sun. 
Blue squares are the total solar gamma-ray flux (solar disk 
-I- IC) within 1.5° of the Sun with only statistical error bars. 
Black dots are the solar-disk-only component from Fermi2011. 
The grey band shows the solar-disk-only component found in 
this work. Green circles are the estimated diffuse background 
within 1.5° of the Sun. The differential point-source sensitiv¬ 
ities for HAWC [72] and LHAASO [73, 74] are shown. 


Utilizing the improved photon statistics, we extend the 
observations to 100 GeV. As in Fermi2011, we find that 
the gamma-ray flux is higher than the central value of the 
SSG1991 prediction by about one order of magnitude in 
1-10 GeV, modulo time variation. In addition, we detect 
the solar-disk component in 10-30 GeV at > 5 cr, and in 
30-100 GeV at > 2 cr. This is the first time the Sun is 
detected above 10 GeV in gamma rays. There are no 
theoretical predictions for solar-disk gamma rays in this 
energy range. As a result, our observations demand fur¬ 
ther theoretical investigation. 

Importantly, we find a significant time variation in the 
solar-disk gamma-ray flux over the analysis period, which 
apparently anticorrelates with solar activity. This is the 
first clear observation of such a time variation, though it 
was hinted at in earlier studies [18, 19]. This variation 
was not theoretically predicted, and its large amplitude 
deserve further investigation. Nonetheless, the anticor¬ 
relation with solar activity indicates that the bulk of the 
solar-disk gamma rays can be explained by cosmic-ray 
interactions in the solar atmosphere and the gamma-ray 
production process is strongly affected by the solar mag¬ 
netic fields. 

Future observations with Fermi and other instruments 
may provide even more information about gamma rays 
from the Sun. For example, the anticorrelation of the 
solar-disk gamma-ray flux with solar activity can be fur- 
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ther confirmed with near-future Fermi data. In addition, 
our robust detection (> 5 a) of the total solar gamma- 
ray flux shows that the Sun is a new and promising 
source for large water-Cherenkov gamma-ray telescopes, 
such as HAWC and LHAASO. Observations from water- 
Cherenkov telescopes can provide important insights on 
the gamma-ray production processes in the TeV range. 

This work lays the observational foundation for our fu¬ 
ture theoretical work, where we will investigate in detail 
how cosmic rays interact with the Sun under the influ¬ 
ence of solar magnetic fields. We will study the multi¬ 
messenger signatures from these high energy processes, 
their implications for solar physics, cosmic-ray physics, 
and new physics. Gamma-ray studies of the Sun are still 
in their infancy, but have already yielded interesting re¬ 
sults. Future observations and the accompanying theo¬ 
retical investigations may uncover even greater surprises. 
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